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ABSTRACT 

The aim of this work is to perform a systematic study of the measures of the 
mass and concentration estimated by fitting the convergence profile of a large sample 
of mock galaxy cluster size lenses, created with the publicly available code MOKA. 
We found that the main contribution to the bias in mass and in concentration is 
due to the halo triaxiality and second to the presence of substructures within the 
host halo virial radius. We show that knowing the cluster elongation along the line 
of sight helps in correcting the mass bias, but still keeps a small negative bias for 
the concentration. If these mass and concentration biases will characterize the galaxy 
cluster sample of a wide field survey it will be difficult to well recover within one 
sigma the cosmological parameters that mainly influence the c — M relation, using 
as reference a 3D c — M relation measured in cosmological ./V-body simulation. In 
this work we propose how to correct the c — M relation for projection effects and 
for adiabatic contraction and suggest to use these as reference for real observed data. 
Correcting mass and concentration estimates, as we propose, gives a measurement of 
the cosmological parameter within 1 — a confidence contours. 

Key words: galaxies: halos - cosmology: theory - dark matter - methods: analytical 
- gravitational lensing: strong 



> 1 INTRODUCTION 



X 



Galaxy clusters represent a very important cosmic labora- 
tory. Their abundance and structural properties are corre- 
lated with the content, the formation and the evolution his- 
tory of the universe. Recent observational analyses have fa- 
vored the dark energy as dominant component, contributing 



to the late accelerate expansion of the universe (Perlmutter 
et al.|[l998 1999[ ), and the dark matter as the second one, 



made up of some kind of non-baryonic weakly interacting 
elementary particle, left over from the Big Bang, this rep- 
resents the concordance ACDM model. Galaxy clusters, the 
most bound and late forming structures of the universe, pos- 
sess a matter content which is compatible with the cosmic 
one ( Ettori et al.|2009 L About 85% of the total galaxy clus- 
ter mass is in form of dark matter, while 15% is attributed 
to baryons: 75% of which is in form of hot and 7% of cold 
gas, other baryon fractions count for the remaining 18%. For 
this reason, cluster number counts represent a very impor- 
tant probe for the nature of dark energy and matter, as like 
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as the study of their mass and concentration, that we will 
present and discuss in this work. Observing them at different 
wavelengths gives the possibility to study their whole con- 
tent: the hot and cold gas, stars in galaxies and the presence 
of dark matter. 

Numerical simulations of structure formation in the 
ACDM framework predict the halo concentration, defined 
as the ratio between the virial and the scale radii of the 
density profile, to be a monotonic function of the host halo 
mass: smaller haloes, forming at higher redshifts than larger 
ones, tend to possess a larger value of the concentration. Re- 
cently, increasing the force and the mass resolution in nu- 
merical simulation has given the possibility to span a very 
large halo mass range down to the size of dwarf- galaxies. 
Dolag et ab[ ([20041); |Neto et al.| ( [2007] ) ; |Macci6 et al.| ( [2007] 
20081; Gao et al. (20081, interpreting the results of different 



numerical simulations, agree in finding a power-law relation 
between the average values of mass and concentration. On 
the other hand, more theoretical interpretations have been 
given to the average c — M relation by Zhao et al. (20091; 
|Giocoli et al.| ( |2012b[ ) , that relate the concentration to the 
time at which the main halo progenitor assembles a certain 
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fraction of their mass. Since observationally the halo forma- 
tion time, tf — defined as the time at which the main halo 



shift z — 0.2 from the Millennium Simulation (Springel et al 



progenitor assembles half of its mass (Lacey & Cole 1993 
Giocoli et al.|2007 \ - is not a direct measurable quantity, we 



need a good and unbiased estimate of the host halo mass and 
concentration to translate them in tf ( Giocoli et al.|2012b |. 
The Cluster Lensing And Supernova survey with Hubble 
(CLASH) project (?), together with X-ray and ground-based 
optical observations, is providing the 2D mass distribution 
of a sample of 25 X-ray selected galaxy clusters. The clusters 
are selected as dynamically relaxed and span a range of red- 
shifts going from z = 0.18 up to z = 0.9. The cluster mass 
density distribution is reconstructed combining both strong 
and weak lensing analyses. The strong lensing analysis is 
performed thanks to the identification of multiple images 
( Zitrin et al. 2009 I allowed by the excellent HST angular 



resolution. The reconstruction of the mass distribution is 
performed assuming that the galaxy cluster light traces the 
matter distribution made up of three components: the clus- 
ter galaxies, a dark matter halo and an external shear to ac- 
count for additional ellipticity in the mass distribution. The 
weak lensing analysis, which probes the mass distribution 
outside the strong lensing region, is conducted using wider 



2009 



ground-based images obtained from Subaru (Umetsu et al. 
~ Okabe et al~1|2010a[ ), combined together with KPNO 



Mayall 4-m imaging and spectroscopy from MMT/Hctospec 
In this case, selected background galaxies are identified from 
the wider cluster images and used to measure the observed 
shear and magnification. Deriving the mass profile simul- 
taneously combining the strong and weak lensing analyses 
allows a very good determination of both the halo mass and 
concentration (Meneghe tti et ah]|2010b" Rasia et al.||2012 |. 
Unlike most other methods that are used to measure the 
mass and concentration of galaxy clusters, like X-ray anal- 
yses that need to assume hydrostatic equilibrium for the 
hot intracluster medium or satellite galaxy velocity distri- 



bution function that needs virial equilibrium (Becker et al 



20071, lensing analyses do not require any assumption of 



the dynamical state of the cluster. The interpretation based 
on the assumption of hydrostatic and virial equilibrium is 
sometimes complicated by the presence of massive and/or 
numerous substructures and by the fact that most clusters 
are not really relaxed. 

Most of the analyses performed on combined weak 
and strong lensing observations tend to find clusters over- 
concentrated with respect to what is expected in a ACDM 
universe. This bias can be due both to intrinsic projection ef- 
fects - main halo triaxiality ( Meneghetti et al. 2007b I and/or 
presence of substructures along the line of sight - and to 



the presence of a massive background structure (Coe et al 



20121. To be able to recover a good and unbiased estimate 



of the dark mass and concentration is important to better 
constrain cosmological parameters. 



Meneghetti et al. ( 2010b I; Rasia et al. (20121, study- 



ing galaxy clusters extracted from numerical simulations, 
have shown that fitting projected quantities, using a |Navarro| 
et al. ( 1997 1 (hereafter; NFW) function, recovers a mass that 
could bo underestimated down to 10 — 15%. The halo mass 
tends to be more biased generally when the 2D projected 
mass is unrelaxed due to presence of nearby structures. 

Recently, |Bahe et al.| ( |2012| | have produced mock ob- 
servations of a sample of massive clusters extracted at red- 



2005 1, to obtain a projected estimate for both halo mass and 



concentration through simulated weak lensing only analysis. 
They showed that the slope and the normalization of the de- 
rived 2D c — M relation tend to be lower than the underlying 
3D ones by 40 and 15 percent respectively. An underestima- 
tion of the mass measured from weak gravitational lensing 



has been also found by Becker & Kravtsov (20111, using 



cluster-size haloes extracted from a numerical simulation. 
Their analysis includes correlated and uncorrelated LSS plus 
gaussian noise in the tangential shear measurement. 

Light travelling from background sources to the ob- 
server is also deflected by the presence of large scale struc- 
tures (LSS) along the line of sight: the observed lensing sig- 
nal is a combination of the cluster and the cosmic shear sig- 
nal (Hockstra 2003). This combination limits the accuracy 
with which both mass and concentration can be estimated, 



not introducing any particular bias. Hoekstra et al. (20111, 



also using the Millennium Simulation results, derived that 
this contribution adds an error budget which is compara- 
ble to the statistical errors. In order to bypass the effect of 
LSS and also the errors and the biases that could be intro- 
duced on the cluster mass reconstructions by the data analy- 
sis software pipelines, we build a large sample of galaxy clus- 
ter convergence maps using the publicly available software 



M0KA (Giocoli et al. 2012a I, and fit their derived spherical 



averaged profile to recover their mass and concentration. 
Throught this work we use a ACDM reference model, 



according to WMAP7-year results (Komatsu et al. 20111: 
fl m = 0.272 (DM+baryons), ft A = 0.728, w = -1, a s = 
0.809 and h = 0.704 at z = 0. 

The paper is organized as follow. In Section [2] we sum- 
marize the ingredients used to create our cluster size-haloes. 
In Section [3] we quantify the bias and scatter caused by in- 
trinsic halo properties. Section [4] is devoted to study how 
well cosmological parameters can be constrained using the 
c— M relation of objects extracted from a wide field survey. 
Summary and conclusions are discussed in Section [5] 



2 THE METHOD 
2.1 The MOKA Code 



M0KA ( [Giocoli et al. 2012a I is a new and fast code able to 
create realistic gravitational lenses, whose size spans from 
galaxies up to clusters of galaxie^] The algorithm uses the 
recent results obtained from numerical simulations to model 
the structural properties of the lenses. In creating cluster 
2D- maps MOKA is very fast, its utilization is friendly allowing 
the possibility of switching on and off different halo proper- 
ties. The user can also easily change and redefine the char- 
acteristics of the main halo, of subhaloes and of the bright 
central galaxy. 



2.2 Presence of Substructures and of a Bright 
Central Galaxy 

In this work we produced a set of simulations creating a 
sample of different galaxy cluster size lenses at various red- 

1 web page: cgiocoli.wordpress.com/research-interests/moka 
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shifts. We have divided the haloes in nine mass bins from 
10 13,5 up to 10 15 ' 5 M e /h, logarithmically equally spaced 
with dlog(M) = 0.25. In each mass bin, with mass M v i r , we 
have generated a sample of 128 lenses, that differ in concen- 
tration c V i r , randomly sampling a log- normal distribution, 
around the average value for that mass, with a given <7i nc , 
and in subhalo abundance. In what follows we summarize 
the ingredients included in the input file of the algorithm. 

• One of the main result obtained studying virialized 
structures in numerical simulations is that haloes possess 
a relatively flat density profile oc r~ near the center and 
steeper oc r~ 3 outside (Navarro et al. 19971. The scale ra- 



dius r s where the slope changes defines the halo concentra- 
tion Cvir = R v ir/r a - For each halo, the smooth halo compo- 
nent is modelled with a NFW profile, to which is assigned a 
concentration c sm such that: 



PNFW,vir{r\M v i r ,C v i r ) = pNFW,s (f | M a 
Ntoi 

^ Psub(r - ri\m su b,i) , 



+ 



(1) 



where the subscript vir, sm and sub refer to the total halo, 
the smooth component and the substructures respectively, Vi 
indicates the distance of the subhalo from the host halo cen- 
ter and p su b(< r) is the subhalo density profile. The equation 
above is solved numerically by minimizing: 



X (Csm) — \ pNFW,vir(r\M vir , C vir ) 



(2) 



PNFW,s(r\M sm , Csm) + ^ paub (r — Ti \ m su b,i ) 
j=l 

• In the standard scenario of structure formation haloes 
collapse as consequence of gravitational instability, with the 
small systems assembling first, in a denser universe, than 



the large ones ( Wechsler et al. 2002 1. This translates in small 



haloes to be typically more concentrated than the large ones. 



In our model we use the Zhao et al. ( 2009 1 relation to take 



into account this effect, which relates the concentration of 
the halo with the time £o.04 at which the main halo progen- 
itor assembles 4% of its mass through the equation: 

1/8 



c vir (M vlr , z) = 4 i 1 + 



t(z) 



3.75io.o 



(3) 



According to this model we estimate the scatter in concen- 
tration at fixed halo mass using the scatter in to.04 from the 
generalized formation redshift distribution derived by |Gio-| 
|coli et al~| ( |2012b| > . In Figure[l]we show the redshift evolution 
of the scatter in concentration as a function of redshift. From 
the figure, we notice that at low redshift the scatter is con- 



sistent with the results of numerical simulations ( Jing 
Jing fc Suto||2002"l [Wechsler et al.||2002| |Dolag"etaT 



Sheth & Tormen 2004 Neto et al 



2000 



2004 



2007) to be of the order 



of 0.23, while at high redshift <ri nc drops down due to the 
contraction of the typical halo merging time scale. However, 
the decrease of the scatter at high redshift does not seem to 
exist in numerical simulation, as recently noticed by Bhat- 
tacharya et al. (2011 1. The difference can be due to the fact 



that numerical simulations are limited in mass resolution 
going to higher z or that the scatter in c V i r at a fixed mass 
is not totally due to the different assembly halo histories. 
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Figure 1. Scatter in concentration cri n c at a fixed host halo mass 
as a function of z due to the scatter in time when the main halo 
progenitor assembles 4% of its mass. The scatter has been ob- 
tained mass averaging the scatter at fixed halo mass of a sample 
of haloes with mass larger than 3.16 X 1O 13 M /A. 



• The improvement of the computational resources has 
given the possibility to increase mass and force resolution in 
numerical simulations. One of the first results revealed by 
this improvement is that haloes are not completely smooth 
but made up of different clumps. Since merging events are 
not totally efficient, these clumps represent the core of pro- 
genitor haloes accreted during the growth of the main halo 



progenitor (Springel et al. 2001b Gao et al. 2004 Giocoli 
et al.|p008 l. Different studies based on the analysis of nu- 



merical simulations have converged on the fact that the sub- 
halo mass function has a power-law distribution with a slope 
a = —1.9. In modeling the subhalo population in M0KA we 



will use the model proposed by Giocoli et al. (2010a I which 



takes into account the dependence on redshift and on the 
host halo mass of the distribution, plus the scatter at fixed 
mass due to the different assembly history; this means that 
the scatters of the concentration and subhalo mass function, 
for a fixed halo mass, are not independent but are correlated: 
haloes more concentrated, and so forming at higher redshift, 
are more relaxed and typically tend to possess less substruc- 
tures ( |Gao et al.|2004||Giocoli et al. 12008] ). We will populate 
the host halo sampling the subhalo mass function down to 
a mass resolution of 10 10 M@ / h, which is of the order of the 
minimum subhalo mass resolved in current cosmological nu- 



merical simulations (Springel et al. 2005 Gottloeber et al. 
2006| |Boylan-Kolchin et al.|2009[ ) 



The subhalo density distribution inside the hosts tends 
to follow that of the host dark matter component. However, 
the combined effect of gravitational heating and tidal strip- 
ping does not allow subhaloes to be present in the central 
region of the cluster, mainly dominated by the density pro- 
file of the smooth component that increases like r~ . This 
translates into a clump density distribution which is less con- 
centrated and shallower than the NFW profile toward the 
centre. To spatially locate the subhaloes within the virial 
radius of the cluster we will adopt the relation by |Gao et al.| 
(2004), since it takes into account these effects being ex- 
tracted from simulated cluters. 
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• Haloes trace the dark matter density field in the Uni- 
verse and their clustering is a function of both the mass and 
the redshift. Their mutual gravitational interaction tends to 
pull and stretch their matter distribution making the haloes 
triaxial. We model the halo triaxiality using the results by 



Jing & Suto ( 2002 1 from numerical studies. This model cor- 
relates the axial ratios e a /e c and e^/e c (e a < ej < e c ) with 
the typical collapsed mass M* , at the considered redshift, 
and the matter content of the universe f2 m ; apart from this 
no other tighter correlation between the halo triaxiality and 
its formation time is assumed. 

• Orbiting in the potential well of the host, subhalo den- 
sity profile tends to be modified since (i) the particles are 
heated at the local temperature of the host, (ii) pushed to- 
ward the tidal radius and (iii) eventually stripped. This phe- 
nomenology modifies the initial NFW profile into something 
that is well described by a singular isothermal sphere (SIS) 
(Metcalf & Madau 20011. The radius of the subhaloes is 
computed in such away preserve the subhalo mass m su b and 
is defined as: 



Gm s 



2a?, 



(4) 



As discussed in Giocoli et al. (2012a I we notice that typi- 
cally, if the subhalo is located at a distance r from the host 
halo center, the following condition is valid Rt(r) ^ R S ub, 
with Rt(r) defined as ( Tormen et al.|1998 ): 

D _ I 1 

* r \[2- 8 In M mr {r) /d In r]M vlr {r) J ' lj 

where M vir (r) represents the mass density profile of the host 
halo. The SIS profile will be used to model the matter den- 
sity distribution around clumps. 

• In the standard scenario \ White fc Rees| [1978] ), the 
structure formation process is driven by the dark matter 
component. Baryons, feeling the halo gravitational poten- 
tial, shock, cool and eventually form stars. This phenomenol- 
ogy drives the formation of a massive and bright galaxy at 
the center of the main halo progenitor. The correct model- 
ing of the bright central galaxy (BCG) is important since 
it influences the lensing cross section of galaxy clusters 
( Mcncghctti et al. 2003). The central galaxy population 



is modeled in M0KA using the halo occupation distribution 



(HOD) approach. We use the relation by Wang et al. ( 2006 1 
to estimate the BCG stellar mass Mbcg, assuming that 
Mbcg correlates with the depth of the potential well of the 
halo within which it formed, thus: 



Me 



2M„ 



(6) 



{M mr /M )- a + (Mvir/M )-P ' 

In this relation, we include a Gaussian scatter in Mbcg 
for given host halo mass with <jm bcg = 0.148, we set the 
parameters a = 0.39, (3 = 1. 96, log(M a t ar , o) = 10.35 and 
M Q = 3.16 x 10 u M a /ft (see | Wang et al.| ( |2006) for more 
detalis about the parameters). 

• The BCG stellar density profile is described following 



the Hernquist (1990) profile 



pstar(r) — 



Pa 



(r/r f )(l + r/r g 



(7) 



where the scale radius r g is related to the half-mass (or effec- 
tive) radius R e by r g = 0.551# e . As done by |Keeton| ( |2001[ ) 
we define the effective radius to be R e = 0. 003 R v i r - The 



scale density p g can be estimated by the definition of the 
total mass of a Hernquist model: 

M K tar 



Pa 



(8) 



We checked the robustness of our results using a Jaffe ( 1983 \ 
profile which keeps unchanged our results. 

In TableQ]we summarize the whole set of simulations we 
have performed: the first two samples, SPH and SPHwBCG, as- 
sume the haloes to be spherical, while the last two, ELLwBCG 
and ELLwBCGwADC, consider the haloes to be triaxial. In the 
latter case the axial ratios are obtained randomly sampling 



the distributions suggested by Jing & Suto (2002). We also 



notice that the last three samples include the presence of a 
BCG. In all samples we included the contribution from the 
subhalo population. We recall the reader that the total halo 
mass M V i r is defined as the sum of all components, so: 



h + m au b,i [ +Mbcg] ■ 



(9) 



In SPHwBCG and ELLwBCG, the presence of a dissipative bary- 
onic component in the host halo centre does not influence the 
dark matter density distribution. This phenomenology, act- 
ing on the dark matter component mainly in proximity of the 



host halo centre and known as adiabatic contraction (Blu 



menthal et a!T1[T986] |Rix et al.|[l997| |Keeton|[200Tj |Gnedin 
et al.|20TT k The effects of baryons on the total mass profile 
are generally found to only modify halo concentrations at 
the 10% level, although some studies suggest that low mass 
systems may me significantly affected by baryonic cooling 



(Fedeli 20111. ELLwBCGwADC includes adiabatic contraction, 



and in Appendix [^] we discuss how the mass and concen- 
tration estimates are biased with respect to the ELLwBCG 
sample. 

As example, in the Figure [2] we show the spherical av- 
eraged convergence profiles and the maps of two clusters, 
with substructures and BCG. In the top panels, we present 
the case of a spherical system while in the bottom panels 
we show a triaxial one (the label reports the values of the 
axial ratios projected in the plane of the lens and along the 
line of sight). We underline that in all our simulations we 
do not consider perturbations of the projected matter den- 
sity due to the large scale structures along the line of sight 



( Hoekstra 2003 1 . This is because we are mainly interested 



to understand how the halo properties their-selves affect the 
recovered mass and concentration fitting a projected quan- 
tity. The map of each halo is resolved with a resolution of 
1024 x 1024 pixels up to the virial radius R V i r defined as 
the radius enclosing the density contrast predicted by the 
spherical collapse model (Eke et al. 19961. The profile is 



built binning linearly the convergence on annulus of 8 pixels 
on which also the rms is estimated, showed by the shaded 
region in the convergence profile of the figure. 



3 FITTING THE CONVERGENCE PROFILE 

The estimate of the total galaxy cluster mass is important 
because it is by far the most accurate predicted halo prop- 
erty from theory. Ongoing and future surveys require an es- 
timate of the mass with a bias smaller than 10% in order to 
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Table 1. List of Simulation sets performed with M0KA 1.5 



Simulation set name Smooth component Substructures BCG Adiabatic Contraction 



SPH 

SPHwBCG 
ELLwBCG 
ELLwBCGwADC 



spherical 
spherical 
triaxial 
triaxial 



yes 
yes 
yes 
yes 



no 

yes 
yes 
yes 



no 
no 
no 

yes 
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Figure 2. Top panels: spherical averaged convergence profile of a spherical halo with substructures and BCG created with M0KA. The 
shaded region shows the rms of the average convergence in each annulus 8 pixels large. The right subpanel shows the convergence map 
of the central region. Bottom panels: spherical averaged convergence profile of a triaxial halo with substructures and BCG. In the label 
we show the values of the components of the axial ratios in the plane of the lens and along the line of sight. The solid and dashed line 
show the input and the best fit converge profile respectively. The right subpanel shows the convergence map up to half virial radius of 
the triaxial halo. 



allow the usage of cluster counts as cosmological probe. How- 



Okabe et al.| |2010bJ; |Oguri et al.| ([2012]) showed that 
the combination of mass and concentration estimates, from 
weak and weak+strong lensing measurements respectively, 
can be used as an additional cosmological probe through the 
c — M relation. These results have revealed a good agree- 
ment with the theoretical predictions for massive clusters 
but a steeper dependence of the concentration on the mass 
wich what pure iV-body simulations predict. As discussed 



by Fedeli (20111 (see also Appendix |A[) , the steepening of 
the c — M relation reflects the response of the dark matter 
density distribution to the presence of the cold baryons in 
the central regions of haloes, that is stronger in less mas- 
sive, then in more massive haloes, since baryon cooling and 
star formation are more efficient in those systems ( |Conroy| 
et al.|[2007a Andrcon 2010). In this direction it is very im- 



portant to understand the scatter and the bias in the mass 
and in the concentration estimated when fitting a projected 
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quantity: the convergence profile. We are aware of the fact 
that the convergence is not a direct observable and that 
the most standard way to measure the cluster mass profile 
with weak lensing is the tangential shear, but our work is 
meant to study the scatter and the bias introduced in the 
mass and in the concentration measurements going from a 
3D to a 2D quantity of the measurement halo profiles. Since 
our idea is to illustrate the intrinsic scatter and the bias in 
the mass and concentration estimates, obtained when fitting 
projected quantities, we present and discuss our results as- 
suming no observational error in the cluster converge map. 
The situation would have been different if we had used the 
tangential shear profile, where the dominant source of noise 
is the intrinsic ellipticity distribution of background galaxies 
cr e « 0.3 ( jOkabe et al.|2010a| ) that gives a Gaussian noise in 
each region of area A where is measured given by: 



rig a iA 



(10) 



where n ga i represents the background galaxy density dis- 
tribution, that goes from 10 gal/arcmin 2 for Dark En- 
ergy Survey (DES) or similar to 40 gal/arcmin 2 for deep 
ground based observations like Large Synoptic Survey Tele- 
scope (LSST) or space based observation like Euclid. It is 
worth mentioning that for high quality observations with 
large source densities, the effect of shape noise becomes 
comparable or subdominant with respect to the intrinsic 
source scatters: triaxiality, correlated (haloes and subhaloes 
around 3 Mpc/h from the cluster) and uncorrelated large- 
scale structures (systems more distant than 3 Mpc/h from 
the cluster). To isolate the effects of halo substructures and 
triaxiality we will describe our results using the cluster con- 
vergence profile. We also underline that the results are spe- 
cific to the method chosen for measuring the halo properties 
(the converge profile). In order to compare the rms of the 
convergence with the statistical observational noise, we show 
in Figure [3] the convergence profile of three clusters located 
at redshift z — 0.03, 0.23 and 0.47. In each panel the shaded 
region represents a K while the solid curves enclose the statis- 
tical observational noise a 3 . We have considered 35, 30 and 
25 gal/arcmin 2 for the cluster located at z = 0.03, 0.23 and 
0.47, respectively. From the figure we notice that while in the 
inner region of the clusters the two errors are almost identi- 
cal, going toward the virial radius the noise expected from 
the intrinsic ellipticities of source galaxies becomes larger 
than the rms of the convergence. 

For each halo the spherical averaged convergence pro- 
file is best fitted using a single component NFW function in 
order to obtain a simultaneous estimate for both the mass 
Mfn and the concentration Cf it . These are obtained consid- 
ering that, if we define x and y as the coordinates on the 
lens plane, the convergence k(t) of a NFW profile can be 
written as ( |Bartelmann||1996| |Meneghetti et al.||2003[ ): 

2p s r s 



-{x,y) 



where r = \/x 2 + y 2 /r s , and F(r): 



F(r) 



: arctan 



= arctanh 



-F(r) 



r + 1 
r - 1 
r + 1 



(11) 



r > 1, 

r < 1, 
r = 1; 



r s is the scale radius, defining the concentration cw = 
Rvir/r„, and p s is the dark matter density at the scale ra- 
dius, 



R/Iyir 



ln(l + Cvir) — 



1 ~\~~ Cvir 



and so 



KNF\v(r) = 



(12) 



(13) 



where E crit = c 2 D s / '(4ivGDiDi s ) is the critical surface mass 
density, depending on the angular-diameter distances Di, 
D s and Di s from the observer to the lens, to the source, and 
from the lens to the source, respectively. From the equations 
above, since knfw(t) depends both on halo mass and con- 
centration, Mfu and Cfu are obtained by minimizing the x 2 
function: 



N ra 



X \Mf it ,c f u) = 



[K NF w(ri\M fu, Cfit) — 



(14) 

where N ra dial bins represents the number of radial bins on 
which the convergence profile is estimated. In the default 
setting wich considers the whole profile [0, l]R V i r the number 
of radial bins is 64, since the map is made by 1024 x 1024 
pixels and the profile is linearly built on annuli of 8 pixels. 



3.1 Spherical Haloes 

Figure [4] shows the rescaled mass Qm = Mfu /M v t r and 
concentration Qc = Cf it /c vir estimates - where M vir and 
c V i r represent the true halo mass and concentration - as a 
function of M V i r for the samples SPH (open diamonds) and 
SPHwBCG (filled circles) at six different redshifts. The esti- 
mates have been obtained by fitting the profile up to the 
halo virial radius: this in practice would mean that for a 
realistic cluster we would need to combine both strong and 
weak lensing data. The different data points with error bars 
represent the median and the two quartiles, respectively of 
the distribution in each mass bin. From the left figure we 
notice that the halo mass is very well recovered without any 
bias and with an error of only few percents. The small trend 
with the halo mass is due to the presence of substructures 
that contaminate the projected mass distribution: the most 
massive haloes possess a larger fraction of mass in substruc- 
tures, because of their later formation time, with respect to 



the least massive systems ( 


Gao et al.||2004 De Lucia et al. 


2004 van den Bosch et al. 


20051. Unlike the mass, the es- 



timated concentration (right panels) shows a negative bias 
and a decreasing behavior with the host halo mass, and so 
with the subhalo abundance. From this figure we notice that 
the concentration estimate is more influenced by the pres- 
ence of substructures in the 2D map than the mass estimate. 
At the smallest masses, the BCG seems to play a more rel- 
evant role, adding an additional bias of few percent. 

To discuss the importance of the radial range used when 
fitting the convergence profile, in Figure|5]we again show the 
estimated mass and the concentration as a function of the 
true halo mass, for the same two halo samples discussed be- 
fore. In this case we show how the estimated masses and 
concentrations depend on the radial range on which the fit 
is performed. In the upper panels, we fit the profile between 
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Figure 4. Median rescaled mass (left panels) and concentration (right panels) estimates as a function of the halo mass, for spherical 
substructured haloes at six different redshifts. The mass of each halo has been estimated by best fitting the spherical averaged convergence 
profile of each halo up to the virial radius R v i r , with a NFW function varying both the mass and the concentration. The two different 
data points refer to haloes without (open diamonds) and with (filled circles) the bright central galaxy. The error bars enclose the first 
and the third quartiles. 



[RmimRmax] /Rvir = [0,1], [0,0.5], [0,0.1] (upper panels of 
Fig. [HJ, [0.1,1], [0.25,1] and [0.5,1] (lower panels). Results 
are shown only for haloes at redshift z = 0.23, since we 
verified that they do not depend on redshift. It is worth 
noticing that we do not observe significant differences apart 
from the following cases: (i) in the fit over the radial range 
[0.25, l]R V ir, where the scatter is almost double; (ii) when 
the region [0, 0.1]R V i r is considered: it roughly corresponds 
to the strong lensing region where the contribution from the 
BCG is comparable to that of the dark matter component, 
while in presence of a BCG the mass is slightly underes- 
timated almost the opposite is valid for the concentration 
(see the different data points) ; (Hi) in the fit over the radial 



range [0.5, l]R V i r , where both mass and concentration are 
overestimated by a factor of 1.5 and 3.5, respectively. 

In order to quantify how much the presence of substruc- 
tures in host haloes biases the mass and the concentration 
estimates, we have divided our sample in haloes with more 
and less substructures. In each mass bin we estimate the me- 
dian host halo mass fraction in substructures (fsub) and we 
define a halo as "substructured" if its fsub is larger than the 
median, f su b > (f su b)- On the contrary, "unsubstructured" 
haloes have / su & < {fsub)- In Figure [6] we show the esti- 
mated mass and concentration at different redshifts (as in 
Figure |4| for the samples of substructured (open triangles) 
and unsubstructured haloes (filled circles). The haloes are 
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Figure 5. Median rescaled mass (left panels) and concentration (right panels) estimates as a function of the halo mass, for spherical and 
substructured haloes at redshift z = 0.23. The top panels show the relation obtained fitting the spherical averaged convergence profile 
up to R v i r , 0.5R v i r and 0.1R v i r ; while the bottom ones from R v i r down to 0.1R v i r , 0.25R„i r and 0.5R v i r . Symbols and error bars are 
as in Figure [4] 




Figure 6. Median rescaled mass (left panels) and concentration (right panels) estimates as a function of the halo mass for the SPHwBCG 
subsamples with more (open triangles) and less (filled circles) substructures. The error bars correspond the first and the third quartile. 



assumed to be spherical and with a BCG at their centres. 
From the figure we notice that, while the estimated masses 
are larger than the true ones for the substructured sample 
than for the unsubstructured, the opposite is valid for the 
estimated concentrations. 



3.2 Triaxial Haloes 

The spherical approximation used to describe the halo mat- 
ter density distribution is far from being realistic. Differ- 
ent studies based on numerical simulations showed that the 



halo shape is triaxial (Sheth et al. II2001 Jing fc Suto||2002 
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Kazantzidis et al.|2004 >. This is because during the collapse 
phase a halo is influenced by its surrounding tidal field which 
stretches and shears its matter distribution. Various obser- 
vations of galaxy clusters, both in optical and in X-rays, 
have confirmed this picture and have revealed that a dis- 
crete number of X-ray-selected clusters present also their 
major axis elongated along the line of sight. This condition 
is typically sufficient, but not necessary, to cause, in their 
optical images, the presence of very distorted gravitational 
arcs and large Einstein rings, and also a discrepancy between 
galaxy cluster masses and concentrations determined from 
X-ray and gravitational lensing observations. Combining X- 
ray, weak and strong lensing observations, Morandi ct al. 
( 2010 1 have proposed a method to determine the galaxy clus- 



ter elongation along the line of sight (with an error of the 
order of 5%) and reconcile the discrepancy of the mass and 
the concentration estimates. The observational estimate of 
the elongation of the cluster halo along of the line of sight 
is very challenging and only possibile for clusters that are 
very massive and located at low redshifts. 

Let us now consider the case of a triaxial prolate or 
oblate halo, and define its and e c ; if the halo 

is spherical in the plane of the sky (e a = ei = 1) and e c is 
parallel of the line of sight, the convergence of the halo will 
be defined as: 



J —! 



p(x, y, C') d C • 



(15) 



with £' = C/ e c- 

In MORA we have included the model proposed by |Jing El\ 



Suto ( 2002 1 to describe the main halo ellipticity, according 



to which the stellar component of the BCG and the subhalo 
spatial distribution are also perturbed. Using the axial ratio 
probability distributions we estimate e a , et and e c , requiring 
that e a e^e c = 1. Once the axial ratios are known, we ran- 
domly orient the halo choosing a point on a sphere identified 
by its azimuthal and elevation angles; according to this we 
rotate the halo ellipsoid and deform the convergence map. 
The elongation e z represents the largest component of the 
ellipsoid axes projected along the line of sight. 

In Figure [7J we show Qm and Q c as a function of the 
host halo mass, obtained by fitting the convergence pro- 
files of the ELLwBCG sample up to R v i r and at six differ- 
ent redshifts. In fitting the convergence profile of the triax- 
ial clusters, and in estimating mass and concentration, we 
proceed in two ways: (i) in the first case we make use di- 
rectly of kmfw(i~) fro m equations 11 and 13 and minimize 



X (Mf ir ,Cfit) (eq. 14 1; («) in the second case we multiply 



knfw{t) by e 2 , the elongation of the halo along the line of 
sight (assumed without error), and estimate Mfu and Cfu 
again minimizing \ 2 ■ We refer to this last case as mass and 
concentration estimates corrected for elongation. 

In the left pa nels of Fig. [7] we notice that the mass 



bias, found also by Meneghetti et al. (2011 \ and Rasia et al. 
( |2012[ > - who fitted the tangential shear profile using sim- 
ulated observations of galaxy clusters - almost disappears 
when we correct for elongation. In the right panels we show 
the results for the concentration. When we do not include 
the elongation in the calculations, both Qm and Q c are un- 
derestimated of 10 — 20%. Conversely, when we correct for 
elongation, the mass estimate is unbiased, while the result- 
ing concentration is still slightly underestimated, since Qm 



and Q c are not independent of each other, considering that 

1/3 

Cvir ~ M vir they do not scale in the same way but as: 
AQ C = (AQ M ) 1/3 . 

To discuss the importance of the range used to fit the 
convergence profile, in Figure|8]we show the median rescaled 
mass and concentration as a function of the true halo mass, 
when the fit of the convergence profile is done on different 
radial ranges. These are the same chosen for the analysis 
presented in Fig. [5] Again we note that the error bars be- 
come larger when fitting on a narrower radial range. In the 
figure we do not observe significant differences in the radial 
ranges [0,1]^, [0,0.5]^, [0.1,l]R vir and [0.25, l]iW 
Fitting the most central region, where the BCG dominates 
the convergence profile, the masses and the concentrations 
are slightly underestimated and overestimated, respectively. 
As show for the spherical case, fitting the profile outside 
0.5R v ir we overestimate both the halo mass and concen- 
tration with not much difference between the two cases in 
which we do or do not correct the convergence profile for 
elongation. 

Fitting the profile without any assumption on the value 
of the elongation creates a strong bias in the mass estimates 
due to the elongation along the line of sight. This can be bet- 
ter quantified dividing the halo sample, in each mass bin ac- 
cording to the level of elongation, in more and less elongated 
systems. In Figure [9] we again show the median masses and 
concentrations as a function of the true mass for the samples 
of haloes with large and small elongation. The separation of 
haloes in these two classes was made by measuring the me- 
dian elongation (e z ) in each mass bin. We define as "less 
elongated" haloes those with e z < (e z ), and vice- versa. This 
analysis shows that the two samples have a large differences 
both in mass (left panel) and in concentration (right panel) 
of about 35 — 40%. We can summarize that while more and 
less substructured haloes have a difference in the estimates 
of the order of few percent, the estimates between more and 
less elongated haloes can differ up to more than 50%, in the 
case of very massive clusters. 

This can originate very important systematics when 
trying to recover cosmological information from the mass 
and concentration estimates based on lensing analysis with- 



out taking into account possible projection effects (Sereno 
|fc Zitrin|2012| |Coe et al.|2012[ ). 



4 COSMOLOGY FROM THE C M RELATION 
OF A WIDE FIELD SURVEY 

Several studies aimed at constraining cosmological parame- 
ters by measuring the c — M relation of galaxy clusters. For 



example, Ettori et al. (2010 1 used a sample of 44 X-ray lumi- 



nous clusters observed with XMM-Newton. They measured 
the masses and the concentrations by assuming spherically 
symmetric X-ray emitting gas distributions in hydrostatic 
equilibrium with the cluster gravitational potentials, which 
are modeled by NFW functionals. By comparing with the- 
oretical c — M relations, they constrain the degeneracy be- 
tween fl m and Checking the consistency of the observed 
cluster c — M relation with the theoretical expectations in 
the framework of the ACDM cosmology is also one of the 
goals of the ongoing CLASH Multi-Cycle- Treasury program 
of the Hubble Space Telescope (?). Interesting results are 
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Figure 7. Median rescaled mass (left panels) and concentration (right panels) estimates as a function of the halo mass, for triaxial 
substructured haloes with BCG, located at six different redshifts. The mass of each halo has been estimated best fitting the spherical 
averaged convergence profile of each halo, up to the virial radius R v i r , with a NFW function varying both the mass and the concentration. 
The two different data points refer to two cases in which the fit has been performed not knowing (filled diamonds) or knowing (filled 
triangles) the halo elongation along the line of sight. The error bars enclose the first and the third quartiles of the distribution in each 
mass bin. 
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Figure 8. Median rescaled mass (left panels) and concentration (right panels) estimates as a function of the halo mass, for triaxial and 
substructured haloes with BCG located at redshift z = 0.23. The top panels show the relation obtained fitting the profile without any 
assumption on the value of the elongation up to R v i r , 0.5R v i r and 0.1R v i r ; while the bottom ones from R v i r down to 0.1i?„i r , 0.25R v i r 
and 0.5i?„i r . The different data points and error bars are the same as in the Figure[7] 
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Figure 9. Median rescaled mass (left panels) and concentration (right panels) estimates as a function of the hist halo mass, for the 
ELLwBCG sample by fitting the profile without any assumption on the value of the elongation using a NFW-funtion. In each mass bin we 
divided the haloes in more (open triangles) and less (filled circles) elongated along the line of sight. The various panels show the relation 
at six different rcdshifts. 



emerging from these analyses, whose significance is however 
limited by the relatively small sizes of the cluster samples. 

The situation is going to significantly improve in the 
near future. Indeed, several optical surveys are planned, 
which will deliver high quality data suitable for the weak 
lensing analysis of la rge samp l es of g alaxy clusters (e.g. 
KiDS, DE!f| LSST), |Kuijken| ( |2010| ); " 



LSST Science Col- 



laborations et al. (2009 1). In particular, in the wide survey 



that will be operated by the Euclid mission, several dozen of 
thousands of galaxy clusters will be detectable up to redshift 
2. Now, we attempt to quantify the robustness of the cos- 
mological constraints that it will be possible to derive from 
the analysis of such large cluster samples. 

The first step towards our goal is to model a cluster se- 
lection function. In this work, we make the simple assump- 
tion that cluster detections will be not contaminated above 
a constant mass limit of 5 x 10 13 Mg//i for z < 0.4 and 
1.5 x 10 14 M©//i for z > 0.4. This reflects the fact that at 
high redshift we tend to select only massive clusters. More 
detailed selection functions should be derived for specific 
surveys, depending on their depth, observing bands, etc. The 
mass limit used in this work is consistent with the analyt- 
ical estimates for the EUCLID photometric selection func- 
tion ( Laureijs et al.|[2011 1. We use the selection function in 
combination with the Sheth & Tormen ( 1999 1 mass-function 



to generate a population of cluster-sized dark matter halos 
reproducing the mass and redshift distributions expected in 
the context of a corresponding ACDM cosmology (see Sect. 1 
for the exact choice of the cosmological parameters). We use 
these halos to populate a light-cone corresponding to a sur- 



vey area of 15000 sq. degrees. The sampled volume extends 
from z — 0.03 and z — 1.5 and is subdivided in slices equally 
spaced in dlog(l + z). We further use the Zhao et al. ( 2009 1 
c — M relation, including a redshift dependent log-normal 
scatter as shown in Figure [l] to assign a concentration to 
each halo. 

As we discussed in the previous sections, projection ef- 
fects introduce substantial biases in the measurements of 
masses and concentrations. Such biases can be quantified 
using the simulations done with the M0KA code. More specif- 
ically, for any input mass and redshift, we can use the simu- 
lations to derive the expected distribution of the measured 
masses and concentrations, assuming they are obtained by 
fitting the two-dimensional mass profiles of triaxial halos. 
In the left and right panels of Figure |10| we show how the 
rescaled masses and concentrations are expected to change 
as a function of the halo virial mass M vir . The results were 
obtained by combining the ELLwBCG samples at 16 different 
redshifts. Triangles and diamonds refer to the cases where 
the projected mass profiles were fitted with and without 
correcting for the elongation along the line of sight (Eq. 15 
with e c = e z ). They indicate the median rescaled masses 
and concentrations in each mass bin. The error-bars corre- 
spond to the inter-quartile ranges of the distributions in the 
mass bins. We performed linear least-squares fits to the data 
points and we measured the relations given by short- and 
the long-dashed lines in each panel. We can explicit these 
relations as 



log(QAf) 
log(Qc) 



0.025 log(M„ ir ) 
-0.018 log(M„i 



- 0.433; rms = 0.094 (16) 
+ 0.212; rms = 0.081 (17) 



2 http : //www . darkenergysurvey . org 



when no correction for elongation is applied. When the fit 
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Figure 10. Median rescaled mass (left panel) and concentration (right panel) as a function of the true halo mass, combining the results 
of 16 redshift slices between z = 0.03 up to z = 1.5, for the ELLwBCG sample. Diamonds and triangles show the results for the NFW fits 
without and with correction for elongation, respectively. The error bars show the first and the third quartile of the distribution. The long 
and short-dashed lines show the least-squares fits to the data. 



is done correcting for e z , we obtain 

log(Q M ) = 0.013 log(M„ ir ) - 0.190 ; rms = 0.066 (18) 
log(Qc) = -0.023 log(M„ ir ) + 0.308 ; rms = 0.029 (19) 

We use these equations and their r.m.s. to assign to each 
halo in the light-cone values of Qm and Q c . Therefore, we 
compute the "biased" masses and concentrations as 



M fu 

Cftt 



Mtrue X Qm 



(20) 
(21) 



Once we have realistically included the biases caused 
by projection effects, we can proceed to measure the halo 
Cfu — Mfu relation and use it to constrain the cosmological 
parameters. In particular, in this section we focus on as, 
S7 m , and w. 

We show in Figure [Tl] the results obtained in six dif- 
ferent redshift bins ranging between z = 0.03 (upper-left 
panel) and z = 1.5 (bottom-right panel). The solid black 
curve show the input c — M relation by |Zhao et al. (2009). 
As usual, the green triangles and the blue diamonds indicate 
the concentrations and masses recovered with and without 
correcting for elongation. As we can see, the diamonds are 
bised lower then the triangles with respect to the Zhao et al. 
c — M relation. These reflect the results shown in Figure [lOl 
where we see that correcting for e z generally allows to re- 
cover the correct input mass, but it leads to a still slightly 
lower estimate of the cluster concentration. These results 
seem to indicate that a correction for elongation, as given 
in Eq. [15] with e c = e z , should be applied to the cluster 
mass. As explained before when we correct for elongation, 
the mass estimate is unbiased, while the resulting concen- 
tration is corrected but not enough, since Qm and Q c are 
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Figure 11. Concentration-mass relations at six different red- 
shifts. The different points show the average concentrations and 
r.m.s. in different mass bins. The solid curve indicate the input 
c — M relation of Zhao et al. Diamonds and triangles show the 
results obtained by fitting the profiles without and with correc- 
tion for elongation. The short and the long-dashed lines are the 
functionals that best fit the data (see text for more details). 



Cosmology in 2D 13 



not independent of each other, considering that Cvir ~ Af„ 



1/3 
\l/3 



they do not scale in the same way but as: AQ C = (AQm ) 

We fit the triangles and the diamonds with functions of 
the kind 



bias: NFW fit triaxial bias: NFW fit triaxial, knowing e z 
1 I ^L' 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 = 



C2d(M) = c 3D (M) x N M 



(22) 



where c 3 b is the input c — M relation (in our case the model 



by Zhao et al. (2009)) and we find that they are well fitted 



by the following functional forms: 

C 1.630 M~ 018 , (short - dashed) 
C2D{M) = c 3D (M) x I 

{ 2.033 M~ 023 , (long -dashed), 

(23) 

respectively. Finally, we use the recovered Cfu — Mfu rela- 
tions to constrain the cosmological parameters. To do so, 
we compare them to corresponding c — M relations of Zhao 
et al. recalculated by varying the cosmological parameters. 
In this experiment, we vary only two of the investigated 
cosmological parameters simultaneously, while we keep the 
third fixed and equal to value used in our simulations. The 
purpose of this experiment is to quantify the degeneracies 
between cosmological parameters when constraining them 
using the c — M relation. If the couple of parameters which 
is changed is ($, \l/), we define the \ 2 variable as 



redshifts N mass 

E E 



(24) 



{l0g[cf it ,i(Zj)] - \0g[Cmodel(Mfit,i, Zj\$, *)]} 2 



where o"i og c,i represents the concentration rms in a given 
mass bin. Then, we find the couples of cosmological param- 
eters, ($,'5'), that minimize \ 2 ■ 

In Figure [T2| we show the 1,2, and 3<r confidence levels 
for the couples of cosmological parameters (as,w) (top pan- 
els), (w,£l m ) (middle panels), and (as, id) (bottom panels). 
The results refer to a 15000 sq. degrees survey, combining 
measurements of the c—M relation in 16 redshift bins be- 
tween z — 0.03 and z = 1.5. In the left panels , we fit the 
input c—M relation, adding only a log-normal scatter in 
c V i r . As expected, we recover the input sets of cosmolog- 
ical parameters, corresponding to the ACDM cosmological 
model. In the middle column, we use the Cfu — Mfu relation 
measured without applying any correction for elongation. In 
this case, the c — M relation is best fitted by models with 
lower erg and lower Q. m when we fix w — — 1 (upper panel). 
Keeping as — 0.809 (middle panel), the results are sensi- 
tive to Q m , and values of w lower than —1 are preferred 
(phantom dark energy). In both upper and middle panels, 
however, the input cosmological parameters are outside the 
3a of the recovered sets of parameters. Fixing f2 m = 0.272 
(bottom panel), the input values of erg and w are outside 
the measurements. The right panels refer to the case where 
concentrations and masses are recovered by correcting for 
the true halo elongations using Eq. [15] with e c = e z . The 
results of the fit in this case are qualitatively similar to 
those shown in the panels in the middle column, but the 
offset between recovered and input cosmological parameters 
is slightly weaker. 

In Table [2j we quantify the degeneracy between the 
couples of cosmological parameters estimated from the 




Figure 13. 1, 2 and 3a confidence levels in the planes (crg,f2 m 
(top panels), erg — f! m , and erg — w (bottom panels). Here, the 
relations given in Eq. |23| are used to correct the biases on the 
c — M relations.. 



concentration-mass relation. Results are reported for 3 — a 
and 1 — a (in parentheses) confidence limits. 

In Figure [13] we present again the confidence levels for 
the cosmological parameters obtained when fitting the c—M 
relations in the cases without and with correction for elon- 
gation (left and right panels, respectively). In these cases, 
however, we used the relations in Eq. [23] to correct for the 
biases in the c — M relations. As expected the input cosmo- 
logical models are now perfectly recovered within lcr. 



5 SUMMARY AND CONCLUSIONS 

In this work, we studied how well the three-dimensional 
masses and concentrations of galaxy clusters are recovered 
from their projected mass distributions. We did this by 
fitting the convergence profiles of a large sample cluster- 
sized halos generated with the publicly available code M0KA 
adopting NFW-functionals. We investigated several sources 
of possible biases in the measurements. In particular, we 
studied the biases introduced by the presence of 1) sub- 
structures of different masses, 2) a massive galaxy at the 
center of clusters, and 3) cluster triaxiality. Additionally, we 
considered the systematics related to the choice of the ra- 
dial range where the projected mass profile is fitted. We also 
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no bias: only log-normal scatter in c„; r bias: NFW fit triaxial bias: NFW fit triaxial, knowing e z 
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Figure 12. 1, 2 and 3 a confidence levels in the planes crg,f2 m (top panels) and crg,w (bottom panels). The results refer to a 15000 sq. 
degrees survey, combining measurements of the c — M relation in 16 redshift bins between z = 0.03 and z = 1.5. In the left panels, we 
fit the input c — M relation, adding only a log- normal scatter in c v i r . In the middle column, we use the Cfn — Men relation measured 
without applying any correction for elongation. Finally, the right panels refer to the case where, concentrations and masses are recovered 
by correcting for the true halo elongations using Eq. |15| with e c = e z . The true sets of cosmological parameters are indicated by the 
black crosses (ACDM cosmological model). 



Table 2. Degeneracy between cosmological parameters estimated from the concentration-mass relation. The subscript (0) indicates the 
input parameters, which correspond to a ACDM cosmological model (f! m .o = 0.272, wq = — 1 and <T8,o = 0.809) 



parameters 


a 


P 


o"8 _ p ( \ Q 

°"8,0 V^m,0/ 


-0.235 (-0.265) 


0.983 (0.994) 


w _ f n m \ a 


1.082 (1.112) 


-1.019 (-1.013) 




0.231 (0.231) 


0.999 (1.000) 
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investigated how the mass and concentration measurements 
change by adding a correction factor that takes into account 
the cluster elongation along the line-of-sight. Finally, we in- 
vestigated the robustness of cosmological constraints derived 
from the concentration-mass relation measured in large sam- 
ples of galaxy clusters. 

Our findings can be summarized as follows: 

• by distributing substructures in spherically symmetric 
cluster halos, we found that the concentrations recovered 
from cluster convergence profiles are on average biased low 
by ~ 5 — 6%. This bias tends to increase as a function of 
mass, because the most massive haloes are expected to pos- 
sess a larger fraction of their mass in substructures. These 
substructures have an even smaller effect on the cluster mass 
estimates, which differ from the input masses by just a few 
percent; 

• the presence of massive galaxies at the cluster center 
generally causes tiny over-estimate of the mass, when the 
convergence profile is fitted with simple NFW functionals. 
Such over-estimate is slightly more evident for small-mass 
halos, where it is of order few percent; 

• the above mentioned results have little dependence on 
redshift, while they depend more significantly on the radial 
range chosen to perform the fit of the convergence profile. 
More precisely, the effects of massive central galaxies are 
remarkable if the fit is restricted to the very inner region of 
the clusters (r < 0.1R vir ). In this case, the cluster mass can 
be under-estimate by up to 10%, while the concentration 
can be over-estimated by the same amount; 

modeling the halo triaxiality following the Jing fc Suto| 



( 2002 1 prescription and assuming random halo orientations, 



the three-dimensional masses recovered from the fit of the 
convergence profiles are on average biased low by ~ 10—20%. 
However, for more massive clusters the mass bias is not 



larger than 10%, as also found by Meneghetti et al. (20111; 



Rasia et al. ( 2012 I using the tangential shear profile of sim- 
ulated galaxy clusters. This bias is larger than what found 



by Corless & King ( 2009 I whose triaxial clusters are con- 



structed without substructures and with a different model 



for the axial ratios e a /e c and et/e c (Shaw et al. 20061 im- 
posing e a ei,e c = 1, as we do, but with e a < et < e c = 1. We 
explained such under-estimate of the mass as due to a larger 
probability to observe clusters elongated on the plane of the 
sky than along the line-of-sight and to the combined effects 
of triaxiality and presence of substructures - as quantified 
and discussed in Section [3] We notice that in this case, we 
also measure low concentrations, this bias in halo proper- 
ties estimates has also been studied and discussed is several 
other studies ( Clowe et al.||2004 Oguri et al.||2005 Gavazzi 



2005 Corless & King 2007 



• introducing a correction for elongation multipling the 
convergence by e z , we found that the bias in mass can be 
removed. This correction also moves Q c toward unity but 
still keeping a small negative bias; 

• by simulating a cluster survey covering an area of 15000 
sq. degrees, we estimated that triaxiality and substructure 
induced biases on the 3D c—M relation, that would translate 
into biases on the values of cosmological parameters like ag , 
Q m , and w, if no correction for elongation is applied. The 
biases decrease when correcting for elongation. 

We conclude that, in order to use the cluster c — M 



relation, measured from the cluster projected density profiles 
as they might be determined by means of the cluster lensing 
signal (see e.g. Coe et al. 2012 Sereno & Zitrin 20121, to 



constrain the cosmological parameters, it is necessary to take 
into account possible biases due to substructures and halo 
triaxiality. This can be done introducing a 2D c — M relation 
correcting the 3D one for projection effects. 
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APPENDIX A: ADIABATIC CONTRACTION 

The presence of a dissipative baryonic component influences 
the dark matter distribution. This is because cold baryons 
that settle at the host halo centre act contracting the dark 
matter. 



Blumenthal et al. ( 1986 1 described the adiabatic 



contraction analytically, finding a good agreement with nu- 
merical simulations. However the halo contraction, as mea- 
sured in numerical simulations with gas physics, may depend 
on the parameters that model the star formation processes 
and the gas physics. In these cases, the evolution of the 
central density may also be affected by spurious numerical 
effects due to two-body scattering of massive particles. Also 
the differences in some models can be due to a higher density 
threshold for star formation, for example in the case in which 
the supernova energy release is more concentrated and cre- 
ates rapid potential fluctuations near the center. The model 



by Blumenthal et al. ( 1986 I assumes that a spherically sym- 



metric halo can be thought of as a sequence of concentric 
shells, made of particles on circular orbits, which homoge- 
neously contract while conserving the angular momentum. 
The initial and final density profiles - characterized by an 
initial radius ri and a final radius r / , when a central galaxy 
is present - are related by 



where 



r [M B ca(r) + M DM j] = nM D M,i{ r i 



MdMJ = MoM,i (1 — fcool) 



(Al) 



(A2) 



and fcooi is the baryon fraction in the halo that cools to 
form the central galaxy. To solve the adiabatic-contraction 



equation, we need to derive r from equation (All. When 
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the stellar density distribution is described by a Hernquist 
model we can read: 

fcooir 3 + (r + r h f [(1 - f coo i)r - n] m^n) = , (A3) 

while for a Jaffe model we have: 



[fcooi + mi(n)(l - fcooi)} r + 

mi{ri){rj - f c rj + n)r + nrjmi(n) = 



(A4) 



where Vh and r,- represent respectively the scale radius of the 
two models. The two equations have a single positive real 
root. We recall that rrn in the previous equations defines the 
initial mass profile normalized by the halo virial mass. From 
the equation above results that the presence of a dissipative 
baryonic component contracts the dark matter distribution. 

The model we described above, and used to model the 
dark matter contraction in our haloes, is often defined as 
standard adiabatic contraction (SAC). [Gnedin et al.| ( [20rT| 
have defined a modified adiabatic contraction (MAC) model 
that takes into account the results from numerical simula- 
tions in which gas physics and merging events are included. 
They argue that the MAC model gives excellent predictions 
to the data, while the SAC model tends to overestimate the 
dark matter content toward the halo centre. We also stress 
that the adiabatic contraction observed in numerical simu- 
lations depends directly on the adopted parameters. Cosmo- 
logical simulations performed by different authors, with very 
different codes and physics inputs, agree that the contrac- 
tion effect is present, but at a weaker level than suggested by 
the SAC model and with a significant variation from system 
to system and from input model to input model. 

Figure |A1| shows the median ratio of the best fit es- 
timate of the mass (crosses) and the concentration (filled 
circles) between the run with (ELLwBCGwADC) and without 
adiabatic (ELLwBCGwADC) contraction as a function of the 
host halo mass. We notice that in systems where the ratio 
between the BCG mass and the virial mass is higher (small 
systems) the contraction tends to increase the concentra- 
tion up to a 10%, while for haloes where Mbcg / 'M V i r is 
lower (larger systems) the ratio approaches to unity. From 
the figure we also notice that the ratio between the two mass 
estimates is influenced by the adiabatic contraction by only 
few percents. 

In the Figure, the two solid curves, that consider the 
case of a BCG following a |Hernquist| |l990"| ) profile, corre- 
spond 



log 



lo» 



Mfit,wADC 
Mf it} nADC 
Cfit,wADC 



exp [— (log M vir 
-- exp [— (log M vir 



13.58)] /88.80, (A5) 
- 14.12)] /51.04,(A6) 



, Cfit }7l ADC , 

respectively, while the dashed curves, that consider the case 



of a BCG following a Jaffe ( 1983) model correspond 



log 



lo» 



Mfit tW ADC 
Mfit,„ADC 
Cfit,wADC 



exp[ 



exp I 



(log M vir 

- (log M vir 



13.95)] /89.72, (A7) 
- 14.53)] /49.04,(A8) 



. Cfit } nADC . 

which best fit the data points. These equations can be used 
to modify the concentration-mass relation measured in pure 
dark matter simulation to include the effects of a dissipative 
stellar component at the host halo centre. 



1 .2 



1.1 




Figure Al. Median ratio between the concentration and mass 
estimates obtained from the simulations performed at z = 0.03, 
0.23, 0.47, 0.75, 1.09 and 1.50 with and without the inclusion of 
the adiabatic contraction, as a function of the host halo mass. 
We show the results for the BCG following a [Hernquist | |1990| 
and a |Jaffe| | |1983[ l profile. The error bars show the first and the 
third quartiles. The curves show the best fit to the data — see 
equations j|A6| (| A5| \A&\ \A7\. 
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